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Abstract: We present a detailed analysis together with numerical calculations on one-loop 
contributions to the neutron electric dipole moment from supersymmetry without R parity, 
focusing on the quark-scalar loop contributions. Complete formulas are given for the var- 
ious contributions through the quark dipole operators. Analytical expressions illustrating 
the explicit role of the R-parity violating parameters are given following perturbative di- 
agonalization of mass-squared matrices for the scalars. Dominant contributions come from 
the combinations B*\' i j 1 for which we obtain robust bounds. Even if the R-parity violating 
couplings are real, CKM phase does induce RPV contribution and for some cases such a 
contribution is as strong as contribution from phases in the R-parity violating couplings. 
Hence, we have bounds directly on {B^X^-A even if the parameters are all real. 
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1. Introduction 

Forty years after the discovery of CP-violation its experimentally observed effects in 
the K and B-meson systems @ are generally compatible with the standard model (SM) 
predictions with the Kobayashi-Maskawa (KM) phase as its sole source. Electric dipole 
moments (EDMs) of elementary particles are likely to be the best candidate for the evidence 
of CP violating physics beyond the SM. The presence of an EDM itself violates CP@. 
Neutron being electrically neutral, is studied extensively by experimentalists, for a possible 
nonvanshing EDM. In the SM, the KM phase generated EDM of neutron is of the order 
1CP 32 ecm Q, being still way below the current experimental bound of 6.3 x 10 -26 ecm M. 

Many extensions of the SM are expected to give potentially large EDM contributions. 
For instance, in the minimal supersymmetric standard model (MSSM), which is no doubt 
the most popular candidate theory for physics beyond SM, the plausible extra EDM con- 
tribution of substantial magnitude has led to the SUSY-CP problemQ. If one simply 
takes the minimal supersymmetric spectrum of the SM and imposes nothing more than 
the gauge symmetries while admitting soft SUSY breaking, the generic supersymmetric 
standard model (GSSM)Q would result. When the large number of baryon or lepton num- 
ber violating terms are removed by imposing an ad hoc discrete symmetry called R-parity, 
one obtains the MSSM Lagrangian. GSSM is a complete theory of supersymmetry (SUSY) 
without R-parity, where all kinds of R-parity violating (RPV) terms are admitted without 
bias. It is generally better motivated than ad hoc versions of RPV theories. The MSSM 
itself, without extension such as adding SM singlet superfields and admitting violation 
of lepton number, cannot accommodate neutrino mass mixings and hence oscillations (8|. 
Given SUSY, the GSSM is actually conceptually the simplest framework to accommodate 
the latter. The large number of a priori arbitrary R-parity violating terms do make phe- 
nomenology complicated. However, the origin of the (pattern of) values for the couplings 
may be considered to be on the same footing as that of the SM Yukawa couplings. For 
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example, it has been shown in || that one can understand the origin, pattern and magni- 
tude of all the R-parity violating terms as a result of a spontaneously broken anomalous 
Abelian family symmetry model. 

The simplest approach to study quantum corrections to the neutron EDM is to analyze 
loop contributions to the u and d quark dipole moments and obtain the overall neutron 
number based on the valence quark model[10|. The neutron EDM is given by 



d n = -(4d d -du) J] , (1.1) 

where r] ~ 1.53 is a QCD correction factor from renormalization group evolution [pT], 12]. 
In the case that only trilinear RPV couplings are admitted, it has been shown that nonva- 
nishing results come in only starting at two- loop level [13]. Striking one- loop contributions 
are, however, identified and discussed based on the GSSM framework [[IJ, [l5|. Similarly, 
flavor off-diagonal dipole moment contributing to the case of the b —* in |l(| decay 



and that of the [i — > e^f [17] decay at one-loop level have also been presented. That is 
the approach taken here. In Ref . fl4| , the one loop EDM contributions from the gluino 
loop, chargino-like loop, neutralino-like loop are studied numerically in some details. It 
shows that the RPV parameter combination fi*X in dominates with small sensitivity to the 
value of tan/3. The experimental bound on neutron EDM is used to constrain the model 



parameter space, especially the RPV part. The deficiency of Ref. [14], missing details of 



quark(-scalar loop), is made up here. We give complete 1-loop formulas for the type of 
contributions to the EDMs of the up- and down-sector quarks, with full incorporation of 
family mixings. We present numerical analysis of quark-scalar loop contributions from all 
possible combinations of RPV parameters. Besides the familiar fj,*X in , there are a list of 
combinations of the type B*X- k which are particularly interesting and will be the focus 
of this paper. Such a calculation has not been reported before. Hence, we obtain bounds 
on combinations of RPV couplings which are otherwise unavailable. Our study also guides 
plausible experimental searches for RPV SUSY signals at the next generation neutron EDM 
experiments. 

We skip the details of the MSSM contribution to neutron EDM since it has already 



been extensively discussed in the literature [18 
2. Formulation and Notation 

We summarize the model here while setting the notation. Details of the formulation 
adopted is elaborated in Ref.jTj]. The most general renormalizable superpotential for the 
supersymmetric SM (without R-parity) can be written as 



W=e ab 
1 



+ - \ aPk LiL b f M] + \ x> jk u°b°b° k , (2.i) 

where (a, b) are SU(2) indices, (i,j,k) are the usual family (flavor) indices, and (a, (3) 
are extended flavor index going from to 3. In the limit where Xijk, X^ k , X'^ k and \i\ all 
vanish, one recovers the expression for the R-parity preserving case, with Lq identified as 
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Hd- Without R-parity imposed, the latter is not a priori distinguishable from the Li's. 
Note that A is antisymmetric in the first two indices, as required by the SU(2) product 
rules, as shown explicitly here with e 12 = — e 2 i = 1- Similarly, A" is antisymmetric in the 
last two indices, from SU(3) C . 

The large number of new parameters involved, however, makes the theory difficult 
to analyze. An optimal parametrization, called the single- VEV parametrization (SVP) 
has been advocated [jlS|1 to make the the task manageable. Here, the choice of an op- 
timal parametrization mainly concerns the 4 L a flavors. Under the SVP, flavor bases 
are chosen such that : 1) among the L Q 's, only Lq, bears a VEV, i.e. {L t ) = 0; 2) 



h jk(= -W) = ^diag{mi, m 2 ,m 3 }; 3) h%(= X ojk = -\ j0k ) = ^diag{m d , m„, m b }; 4) h u lk = 



j 

s/2 



^CKM 



diag{m u , m c , m t }, where v = y/2 (Lq) and v u = y/2(H u ). The big advantage of 
here is that the (tree-level) mass matrices for all the fermions do not involve any of the 
trilinear RPV couplings, though the approach makes no assumption on any RPV coupling 
including even those from soft SUSY breaking; and all the parameters used are uniquely 
defined, with the exception of some possibly removable phases. 

The soft SUSY breaking part of the Lagrangian in GSSM can be written as follows : 



ajb fc 



+ ^6 



+ Q ] rhl Q + U^rh^ U + D^rhl D + Dm?L + &m\ E + rhl u \H u \ 



+ h.c. 



M - - M , - Mo 
+ ^BB + t±WW + -fgg + h.c. 



(2.2) 



where we have separated the R-parity conserving ones from the RPV ones (Hd = Lq) for 
the vl-terms. Note that L^m|Z, unlike the other soft mass terms, is given by a 4 x 4 matrix. 
Explicitly, wf oo is m^ d of the MSSM case while fh£ 's give RPV mass mixings. 

Details of the tree-level mass matrices for all fermions and scalars are summarized in 
Ref.0. For the analytical appreciation of many of the results, approximate expressions of 
all the RPV mass mixings are very useful. The expressions are available from perturbative 
diagonalization of the mass matrices 0. 



3. The Quark Loop Contributions 

We perform calculations of the one-loop EDM diagrams using mass eigenstates with their 
effective couplings. The approach frees our numerical results from the mass-insertion ap- 
proximation more commonly adopted in the type of calculations, while analytical discus- 
sions based of the perturbative diagonalization formulae help to trace the major role of 
the RPV couplings, especially those of the bilinear type. The interesting class of one-loop 
contributions largely overlooked by other authors on RPV physics consist of a bilinear- 
trilinear parameter combinations. The bilinear parameters come into play through mass 
mixings induced among the fermions and scalars (slepton and Higgs states), while the tri- 
linear parameter enters a effective coupling vertex. The basic features are the same as 
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those reported in the studies of the various related processes [14, [R 
our recently available reports on b 



Among the latter, 



s + 7 1 16] are particularly noteworthy, for compari- 
son. The d quark dipole plays a more important role over that of the u quark, when RPV 
contributions are involved. The b — * s + j diagram is, of course, nothing other than a flavor 
off-diagonal version of a down-sector quark dipole moment diagram. 

The SM Yukawa couplings are real and flavor diagonal, not to say very small for 
the u and d quarks. However, the RPV analogs are mostly not so strongly constrained 
in magnitudes f20|| and are sources of flavor mixings. A trilinear A' - fe coupling couples a 
quark to another one, of the same or different charge, and a generic scalar, neutral or 
charged accordingly. The coupling together with a SM Yukawa coupling at another vertex 
contributes to the EDMs through some scalar mass eigenstates with RPV mass mixings 
involved, as illustrated in Fig. |]. As the SM Yukawas are flavor diagonal, the charged 
scalar loops contribute by invoking CKM mixings 1 . Note that under the SVP adopted, 
the X a j k parameters have quark flavor indices actually defined in the d-sector quark mass 
eigenstate basis. Analytically, we have the formula 



O cll 



Air sin 2 6^ 



^ ^ ^ > lm(^ 'nmi^ 'nmi 



m n=l 



J 7 

M? 



(Q f -Qf)F 4 \ 



M? 



Qf> F 3 



Mr, 

M? 



(34) 



where, for / = u (i = 1), the coefficients C n ' mi are defined by the interaction Lagrangian, 



C u = g 2 V(d n )<S> 



& i^i+C* 1 + 75 

2 



''nmi ^ 



Vim) +h.c. 



(3.2) 



pL* Vd n Trin -W* _i_ 

Vui rrin *r\l* 
"CKM u lm i 

92 



v CKM u lj+2)m 

92 



(3.3) 



and, for / = d (i = 1), the coefficients are defined by the interaction Lagrangian, 



c l 1 ~ 7s a, l + 7s 

nmi 2 '' ■ 



$>(di) +h.c. 



(3.4) 



pL* Vun T/~wi* f\l 

'-'nmi "CKM U lm 1 

92 

V 

pR* Vdj uni* <pl , jki^rnk*q-\l (o r\ 

^nmi — ~~~~ CKM lyxi ~~~ ^CKM U (j+2)m ' \°-°) 
92 92 

with the Yl'm denoting a sum over (seven) nonzero mass eigenstates of the charged scalar; 
i.e., the unphysical Goldstone mode is dropped from the sum, D l being the diagonalization 
matrix, i.e., D l ^M 2 D l = diag{ M- 1 , m = 1-8} and Mj being 8x8 charged-slepton and 
Higgs mass-squared matrix . Here, % = 1 is the flavor index for the external quark while 

1 For the corresponding flavor off-diagonal transition moment, like b — » s + 7 , scalar loop contributions 
do exist 
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n is for the quark running inside the loop. Note that the unphysical Goldstone mode is 
dropped from the scalar sum because it is rather a part of the gauge loop contribution, 
which obviously is real and does not affect the EDM calculation. 

In general, there also exist the contributions from neutral scalar loop (involving the 
mixing of neutral Higgs with sneutrino in the loop). These contributions lack the top 
Yukawa and top mass effects that enhance the contributions of charged-slepton loop, and 
hence less important. The formula is given as: 

v / <f> m n=l Sm \ lm / 

where, for / = u (i = 1), the coefficients M^mi are defined by the interaction Lagrangian, 
C u = g^{u n )<S>\ct? m ) 



n.rn.i. .-■•» 1 * n.m.i. -» 



+h.c. , (3.7) 



92 V2 



1 



92 a/2 



1 



(3.8) 



and, for / = d (i = 1), the coefficients M^mi are defined by the interaction Lagrangian, 



C d = g^(d n )<S>{4? m ) 



n.m.i, i-v n.m.i. . -\ 



*(di) +h.c. 



(3.9) 



* ti vn t 



92 dm V2 [ 2m 



iV. 



7m\ 



jin 



92 V2 



Oi 

92 



1 



P2m + iV 7m \ 



A' 



jni 



52 



V [j+2)m ~ lT> U+7)m 



V (j+2)m + W (j+7)m 



(3.10) 



with the 2^ m denoting a sum over (10) nonzero mass eigenstates of the neutral scalar; i.e., 
the unphysical Goldstone mode is dropped from the sum, D s being the diagonalization 
matrix for the 10 x 10 mass matrix for neutral scalars (real and symmetric, written in 
terms of scalar and pseudo-scalar parts). Again i = 1 is the flavor index for the external 
quark while n is for the quark running inside the loop. 



4. Results 

Quark EDMs involve violation of CP but not imparity. Thus R-parity violating parameters 
should come in combinations that conserve R-parity. From an inspection of formulae it is 
clear that the contributions from two A' couplings cannot lead to EDM as these violate 
lepton number by two units which has to be compensated by Majorana like mass insertions 
for neutrino or sneutrino propagators. If one is willing to admit more than two A' to be 
non-zero, than in principle one can have contributions to EDM from the Majorana like 
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mass insertions but these would be highly suppressed. With only two RPV couplings, the 
only other possibility is to have a combination of trilinear and a bilinear (//j, B{ or w| Q .) 
RPV couplings in such a way that lepton number is conserved. However not all the three 
bilinears mentioned above are independent as they are related by tadpole relation in the 
single VEV parametrization 



Bi tan/3 



™2 



in 



-t-o 



i + MoMi • 



(4.1) 



Using the above tadpole equation we eliminate rh LQ . in favor of /ij and Bi. Contributions 



from the combination /^*A- jfc , through squark loops, had been extensively studied in |Lj 
in detail. Here, we shall focus on the -B*A^- fc kind of combination as illustrated in Fig.|l[ 
Such a combination can contribute through charged-scalar (charged-slepton, charged Higgs 
mixing) loop or neutral scalar (sneutrino, neutral Higgs mixing) loop. The fermions running 
inside the loops are quarks, instead of the gluon or a colorless fermion as in the case of the 
squark loops. Before we discuss the numerical results it would be worth-while to discuss 
the analytical results which can then be compared with numerical results. 

Let us focus on RPV part of Im^^C^J, in Eq. (|3.l[) for the case of d quark EDM. 
It is given as: 



For the u quark dipole, we have 

MCLiCaW = Im 



/* 



•jkl CKM u (j+2)mj 



(4.2) 



(4.3) 



Interestingly, the entries of the slepton-Higgs diagonalizing matrix involved are the same in 
both terms above. Bilinear RPV terms are hidden inside the slepton-Higgs diagonalization 
matrix elements. To see the explicit dependence on bilinear RPV terms let us look at the 
diagonalizing matrix elements of charged-slepton Higgs mass matrix, ^L^w^im) more 
closely. When summed over the index m, it of course gives zero, owing to unitarity. This is 
possible only for the case of exact mass degeneracy among the scalars when loop functions 
factor out in the sum over m scalars. In fact, the final result involves a double summation 
over the m scalar and n fermion (quark) mass eigenstates. Either the unitarity constraint 
over the former sum or the GIM cancellation over the latter predicts a null result whenever 

f M f> \ ( M f> \ 

the mass dependent loop functions F4 1 j4p- J and F3 1 -^f- J can be factored out of the 

corresponding summation due to mass degeneracy. In reality however, these elements are 
multiplied by non-universal loop functions giving non-zero EDM. Restricting to first order 
in perturbation expansion of the mass-eigenstates, ^(j^m^im ^ s non - zero f° r Tn = 2 and 
m = j + 2, both giving similar dependence on RPV but with opposite sign (m = 1 is the 
unphysical Goldstone state that is dropped from the sum here). For m = 2 one obtains @: 



V 



1* 



(i+2)2 ^12 



B* 

~M1 X ° (1 > 



(4.4) 



Here, denotes the difference in the relevant diagonal entries (generic mass-squared 
parameter of the order of soft mass scale) in charged-slepton Higgs mass matrix. To first 
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order in perturbation expansion there is no contribution to EDM from jXi- If one 
second order in perturbation expansion, one gets a contribution from the term m 
which is given as Q: 



goes to 
= j + 5 



,1 -W* 
(j+2)(j+5) U l(j+5) 



Ml 



X 



(A* — fj, tan [3)rrij 
Ml 



V2M wS m(3{fx k X* kjj ) 



(4.5) 



There are a few important things to be noted here: 

1. Notice that ^ enter at second order in perturbation expansion. Moreover they are 
accompanied by corresponding charged-lepton mass-squared and hence /U*A^ - fe contri- 
butions are always suppressed (later below we will elaborate on this). Also notice that 
in principle trilinear parameter X^jj also contribute through fJ-k^kjj- But they have to 
be present in addition to the trilinear parameter X'- k , thus making it a fourth-order 
effect in perturbation and hence negligible. Thus, we will focus on B*X'^ k effects 
which are interesting. 

2. Even if all RPV parameters are real, CKM phase in conjunction with real RPV 
parameters could still induce EDM. As we will soon see, this could be sizable. 

3. It is clear that d quark EDM receives much larger contribution owing to top- Yukawa 
and proportionality to top mass. There are nine trilinear RPV couplings A^ that 
contribute to d quark EDM. A^ 31 has the largest impact owing to the least CKM 
suppression. 

4. All the 27 trilinear RPV couplings (A^- fc ) contribute to u quark EDM. However, the 
absence of enhancement from the top mass in the loop and the uniform proportion- 
ality to up- Yukawa considerably weakens the type of contribution to u quark EDM. 

5. There is no RPV neutral scalar loop contribution to u quark EDM. However, there 
are one-loop contributions to d quark EDM from the neutral-Higgs sneutrino mixing 
due to the combination of BiX' in with Majorana like mass-insertion in the loop 2 . 
From the EDM formula, it is clear that this would be about the same magnitude 
as the u quark EDM due to charged-Higgs charged slepton mixing and hence highly 
suppressed. 

From the above analytical discussion, we illustrated clearly how the various combina- 
tions of trilinear and bilinear RPV parameters contribute to neutron EDM. We now focus 
on the numerical results. In order to focus on individual contributions we keep a pair of 
RPV coupling (a bilinear and a trilinear) to be non-zero at a time to study its impact. 
We have chosen all sleptons and down-type Higgs to be 100 GeV (up-type Higgs mass and 

2 The Majorana like mass-insertion can be considered a result of the non-zero Bi. It manifest itself in our 
exact mass eigenstate calculation as a mismatch between the corresponding scalar and pseudo-scalar parts 
complex "sneutrino" state which would otherwise have contributions canceling among themselves. With 
the non-zero Bi, the involved diagram is one with two X'm coupling vertices and a internal d quark. Hence, 
the type of contribution is possible only with the single A' coupling. 
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i?o being determined from electroweak symmetry breaking conditions), /x parameter to be 
-300 GeV, A parameter at the value of 300 GeV, and tan/3 = 3 (we will show the impact 
of some parameter variations below). Since the A' couplings are on the same footing as 
the standard Yukawa couplings, they are in general complex. We admit a phase of 7r/4 for 
non-zero A' couplings while still keeping the CKM phase. The phases of the other R-parity 
conserving parameters are switched off. Since it is the complex phase of the B*X'^ k product 
that is importance, we put the phase for Bi to be zero without loss of generality. We do 
not assume any hierarchy in the sleptonic spectrum. Hence, it is immaterial which of the 
bilinear parameter Bi is chosen to be non-zero. We choose B% to be non-zero but all the 
results hold good for B\ and B2 as well. 

In Fig.(|2|) we have plotted the neutron EDM versus the most important combination 
|-B 3 A 331 | normalized by /i 2 . As mentioned above, it can be clearly seen that the d quark 
EDM contribution from the charged-slepton loop (green circles) dominates over the u quark 
EDM (blue triangles). The experimental upper bound is shown with the pink horizontal 
dotted line and the QCD corrected total neutron EDM is shown in red squares. From 
the plot one can extract an upper bound of 4.0 • 10~ 5 on the combination of couplings 
I -^3 ^331 1 /A*o- However this bound corresponds to a fixed value of relative phase (vr/4) in 
the RPV parameter combination) as well as other supersymmetric parameters. Later in 
the table [|, we show the effects of variation in the major parameters. With a similar value 
of input parameters we obtain an upper bound of 2.2 • 10~ 4 for \B 3 X' 32 i\/ fJ-o and 8.4 • 10~ 4 
for |5|A 311 |///q. The hierarchy in these bounds is due to the relevant CKM factors. 

Fig. (|3]) shows the contours of neutron EDM in the plane of magnitudes for the couplings 
A 331 and B3 with the relative phase fixed at 7r/4. The contour for present experimental 
bound is shown in pink dotted line. The plot shows that a sizable region of the parameter 
space are ruled out. Successive contours show smaller values of neutron EDM with the 
smallest being 10~ 27 e cm. 

So far we have kept certain parameters like the phase of the RPV combination, the 
H parameter and the slepton mass spectrum fixed. To get a better understanding of the 
allowed regions in the overall parameter space, let us focus on variations of the parameters 
one at a time. In Fig.(Q) we have plotted neutron EDM versus the slepton mass parameter 
% = rfiE (with I S3 1 = 200 GeV 2 , |A 331 | = .05 and relative phase of 7r/4). As expected the 
neutron EDM falls with the increasing slepton mass. More careful checking reveals that 
the result is sensitive only to one slepton mass parameter, the mass of the third L-handed 
slepton here. It is also easy to understand from our analytical formulas that the dominating 
contributions among the various scalar mass eigenstates for the case of B^X' i j k come from 
the i-th. L-handed slepton and the Higgs. In Fig. (|5j) we have shown the variation of neutron 
EDM with the /j, parameter. Although parameter /i does not directly figure in the EDM 
formula, its influence is felt through the Higgs spectrum. Larger /i leads to heavier Higgs 
spectrum which suppresses the EDM contribution. As our analytical formulas show there 
is no strong dependence on tan/3 (we have checked this numerically as well) and hence we 
have kept tan/3 = 3 fixed in all the plots. 

So far we have discussed only the combination B*X' i31 for the purpose of illustration. In 
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table |, we list the bounds on the nine combinations |-B*A^--|J normalized by (100 GeV) 2 . It 
is interesting to note the difference in bounds for case (a) and (b) for coupling combination 
|-B*A^ 31 | and |-B*A^ 21 |. The inputs for Case (b) are identical to case (a) except for RPV 
phase being zero for case (b). Case (b) thus relies solely on CKM phase. Interestingly 
the bound for |_B*A- 31 | changes very marginally from case (a) to (b) whereas the bound for 
1-^*^21 1 weakens by about an order of magnitude. To understand this difference in behavior 
of |-B*A- 31 | and |-B*A- 21 |, with and without a complex phase in the RPV couplings, we have 
plotted in Fig. ^ the allowed region in the plane of relative phase of B^ A 331 and the |I? 3 A 331 1 
(left) and relative phase of i?|A 321 and the l-B^A^I (right). One can see that the bound 
for |-B 3 A 331 | is about the same for relative phase in I3 3 A 331 of or 7r/4 but the bound for 
|-B 3 A 321 | strengthens by about an order of magnitude as the phase in the RPV coupling 
increases from to 7r/4 suggesting a collaborative effect between the CKM phase and the 
RPV phase. For the case of B*X' iU there is no CKM phase involved. In the table we have 
fixed the sign of fi parameter to be negative. In the Fig. |5| it is seen that absolute value 
of neutron EDM is symmetric with respect to sign of /x parameter and hence positive [i 
should give identical bounds. The variation of bounds in table |l| with changes in parameters 
/x and the slepton and Higgs mass very much follows the pattern found in plots discussed 
earlier. 



Mo 



parameter values 



RPV 

phase 



1^31 1 

(ioo GeV) 2 



Normalized bounds 



(ioo GeV) 2 (ioo GeV) 2 



(a) 


-100 GeV 


100 GeV 


100 GeV 


vr/4 


1.6 


10~ 


-4 


4.8 


io- 


-4 


2.0 • 10" 


-3 


(b) 


-100 GeV 


100 GeV 


100 GeV 





1.8 


10~ 


-4 


4.5 


10" 


-3 


N.A. 




(c) 


-400 GeV 


100 GeV 


100 GeV 


vr/4 


5.0 


10~ 


-4 


3.2 


io- 


-3 


1.1 • 10" 


-2 


(d) 


-800 GeV 


100 GeV 


100 GeV 


vr/4 


1.7 


io- 


-3 


1.0 


io- 


-2 


3.7-10" 


-2 


(e) 


-100 GeV 


400 GeV 


100 GeV 


vr/4 


6.5 


io- 


-1 


5.0 


10" 


-3 


1.8 • 10" 


-2 


(f) 


-100 GeV 


800 GeV 


100 GeV 


vr/4 


2.0 


io- 


-3 


1.7 


10" 


-2 


5.8 • 10" 


-2 


(g) 


-100 GeV 


100 GeV 


300 GeV 


vr/4 


4.9 


io- 


-4 


2.1 


10" 


-3 


7.5 • 10" 


-3 


(h) 


-100 GeV 


100 GeV 


600 GeV 


vr/4 


1.1 


io- 


-3 


6.8 


io- 


-3 


2.4 • 10" 


-2 



Table 1: Here we list the normalized upper bounds for several combinations of bilinear and trilinear 
RPV parameters, with some variation in the input parameters. The bounds essentially depend on 
the values of parameters like fh L , ^ and Tfi Hd (ffl>H u and Bq being determined from EW symmetry 
breaking condition), tan/3 has been kept fixed at 3 as EDM has a very mild dependence on tan/3. 

Before we conclude, we would like to briefly comment on two things. In the table [I], 
we have only mentioned the bounds for the combinations l-B^A^I, whereas earlier in the 
text we did mention that in principle all the twenty-seven A^- fc together with bilinear 
couplings do contribute to the neutron EDM. The couplings other than X' i j 1 contribute to 
u quark EDM and all the contributions are proportional to up- Yukawa (in contrast to the 
presence of a contribution proportional to top- Yukawa for the d quark EDM). Strength 
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of the corresponding contributions is substantially weaker (typically by 5 to 6 orders of 
magnitude) and hence do not lead to meaningful bounds. The other thing is about the 
possible \x\ ■ A^- fe contributions. It can be seen in Eq. (|4.5[ ) that these are second order in 
perturbation and also accompanied by lepton mass /ij. Thus, a typical /^A^ contribution 
is substantially weaker than the corresponding B^X' i j k contribution. For instance, with the 
similar inputs for other SUSY parameters as described above, if one takes fi 3 = 1CT 3 GeV 
(dictated by requirement of sub-ev neutrino masses), Ag 31 = .05 and a relative phase of 
7r/4, one obtains neutron EDM of 6.0- 10~ 32 which is about six orders of magnitude smaller 
than the present experimental bound on neutron EDM. If one goes by the upper bound on 
the mass of v T of 18.2 MeV [22|, could be as large as 7 GeV for tan/3 = 2 and sparticle 
mass ~ 300 GeV @. For (j, 3 = 1 GeV we obtain neutron EDM of 6.1 • 10~ 29 , still about 
three orders of magnitude smaller than the experimental bound. These numbers can be 
easily compared with the At*A^ 11 contributions to the d quark EDM through chargino loop 
in Table 1 of ref.[14|. There the corresponding contribution (with fi 3 = 1 GeV) is much 
weaker (of order 10~ 32 ) as it lacks the top Yukawa and top mass enhancements. In the 
same table one finds that corresponding contribution to gluino loop is much stronger (of 
order 10~ 25 ) due to proportionality to gluino mass and strong coupling constant. In the 
light of above reasons one can appreciate that the contributions due to soft parameter Bi 
are far more dominating in the present scenario of quark loop contributions. 



5. Conclusions 



We have made a systematic study of the influence of bilinear and trilinear RPV couplings 
on the neutron EDM. A combination of bilinear and trilinear RPV is the only way RPV 
parameters can contribute at one loop level. RPV couplings contribute to both, it quark 
as well as d quark EDM. Such a contribution could come from charged-slepton Higgs 
mixing loop (contributing to both u quark as well as d quark EDM) or sneutrino-Higgs 
mixing loop (viable only for combination of nonzero -BjA- n giving Majorana-like scalar 
mass insertion) contributing to only d quark EDM. The class of diagrams all have a quark 
as the fermion running inside the loop. In our analytical expressions obtained based on 
perturbative diagonalization of the scalar mass-squared matrices, we demonstrated that 
charged-slepton Higgs mixing loop contribution to d quark EDM far dominates the other 
contributions due to a diagram with the top-quark in the loop. In our numerical exercise 

\B* -X' | 

we have obtained robust bounds on the combinations - — Vttn; for i, j = 1, 2, 3 that have 

(ioo GeV) 2 ,J ' ' 

not been reported before. Even if the RPV couplings are real, they could still contribute 
to neutron EDM via CKM phase. For some cases CKM phase induced contribution is 
as strong as that due to an explicit complex phase in the RPV couplings. There also 
exist contributions involving /^*A^- fc . However these are higher order effects which are 
further suppressed by proportionality to charged lepton mass. Since [ii are expected to be 
very small (of order 10 -3 GeV) for sub-eV neutrino masses, such contributions are highly 
suppressed. Our results presented here make available a new set of interesting bounds on 
combinations of RPV parameters. 
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Figure 1: d quark EDM due to B*X' i31 combination. Due to top Yukawa and top-mass dependence 
this is the most dominant contribution. 
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Figure 2: (color online) The neutron EDM (in units of e.cm) versus the combination |£? 3 A 331 | 
normalized by The charged scalar contribution (green dotted line) to the d quark EDM is 
clearly seen to dominate over the charged-scalar u quark EDM (blue dotted line). The total QCD 
corrected neutron EDM is shown in red line. The experimental upper bound is shown in pink 
dotted line. 
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Figure 3: (color online) A contour plot for various values of neutron EDM in e.cm in the plane 
of real A 331 and _B 3 . The red line is the present experimental bound for neutron EDM. 
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Figure 4: Neutron EDM vs. the slepton mass parameter ifiL = tue (with |£?3| = 200 GeV 2 , 
1^331 1 = -05 and relative phase of 7r/4). Horizontal line is the experimental bound. 
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Figure 5: Neutron EDM vs. /i parameter with \B 3 \ = 200 GeV 2 , IA33-J = .05 and relative phase 
7r/4. Horizontal line is the experimental bound. 
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Figure 6: Dotted region is the parameter space allowed by neutron EDM upper bound in the plane 
of relative phase of X' 331 B 3 and the | -^331 ^3 1 (left) and relative phase of X' 32 iB 3 and the | A321-B3 1 
(right). 
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